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^ 1 Abstract 



The invention of supersymmetry, almost exactly 25 years ago Q, changed 
the face of high-energy physics. The idea that the observed low-energy gauge 
groups appear due to the process of spontaneous breaking of a single unifying 
group G is also quite popular. The synthesis of these two elements results 
in supersymmetric grand unification. I present (perturbatively) exact results 
regarding the supersymmetric evolution of the gauge couplings from the scale 
of their unification to lower scales. In particular, it is shown how the heavy 
mass thresholds can be properly taken into account to all orders. 



* Extended version 



When I was kindly invited by Prof. Mohapatra, almost a year ago, to give this 
talk, it was supposed that the main subject would be the a s crisis 0. Since then 
the situation with low versus high-energy determinations of a s (M z ) was reviewed 
several times ||, and nothing exciting happened during this year. The irreconcil- 
able contradiction between the values of a s (Mz) from the low-energy measurements 
and the high energy determinations is not going to disappear - this fact gradually 
becomes obvious to everybody. Perhaps, it is worth mentioning a very recent record- 
breaking analysis of deep inelastic scattering [|J which achieves three-loop accuracy, 
for the first time ever in this problem. The analysis confirms, with better stability, 
the earlier conclusion, that a s (Mz) is equal to 0.110 with the uncertainty close to 
0.005. According to Ref. 0], this is a conservative estimate of the error bars. If we 
compare this result with the values of a s extracted from the measurements at the Z 
peak we observe a discrepancy of 2.5 to 3 standard deviations. Thus, the question 
is not whether or not new physics is already with us but, rather, what kind of new 
physics manifests itself in the existing data. Quite a few talks at this Workshop were 
devoted to various new physics scenarios, expressing various degrees of confidence in 
the scenarios under consideration ||. I have nothing to add in this respect. In the 
absence of fresh ideas concerning the a s crisis, I decided to discuss today another 
topic which, on one hand, is rather closely related to the a s problem and, on the 
other hand, did not receive due attention previously - a very peculiar pattern of the 
supersymmetric evolution of the gauge couplings. 



1 Outlining the Problem 

The bold idea of grand unification, after two decades of its existence ||, has become 
an indispensable part of phenomenology. It has obviously attractive features, e.g. 
the fact that it naturally explains the electric charge quantization. Its most unpleas- 
ant feature, to my mind, is the prediction of a great desert above the supersymmetry 
scale (100 GeV range) up to the grand unification scale (10 16 GeV range), an energy 
interval spanning fourteen orders of magnitude or so, where nothing happens. If the 
great desert does indeed exists, high energy physics will practically cease to exist as 
the experimental science we know today, after the discovery of the superpartners. 
I would not like to believe in such a frightening future, and, yet, I am going to 
discuss the supersymmetric grand unification, for two reasons. First, I will empha- 
size some general aspects of supersymmetric evolution of the gauge couplings which 
may well be of importance in other problems, not necessarily related to geography 
of the great desert. Treating the gauge coupling evolution in a standard way 0, 
one overlooks elegant properties which are specific to supersymmetry. My task is to 
reveal these elegant properties. Second, supersymmetric grand unification per se is 
an intellectual game occupying the minds of many, probably, because it predicts a 
relation between the low-energy gauge couplings which is surprisingly close to the 
one observed experimentally. It seems worth joining the club of unifiers merely to 
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understand whether the experimentally observed relation is actually reproduced in 
the most naive and straightforward version of supersymmetric grand unification. 

In the first part of my talk I will focus on the generalities of the supersymmetric 
evolution of the gauge couplings. Assume that there is a (simple) gauge group G 
which is spontaneously broken down to several subgroups. Far above the symmetry 
breaking scale we start from a single gauge coupling. Then it splits, and far below 
the breaking scale we deal with a whole set of independent couplings which, however, 
remember their common origin. The common origin imposes certain constraints. 

There is nothing new in this formulation, of course. Getting observable conse- 
quences requires evolving the constants from the unification scale down to a lower 
scale - this fact was recognized almost immediately || after the emergence of the 
grand unified theories (GUT). Huge amounts of literature are devoted to this is- 
sue. Let me mention early works on supersymmetric grand unification ||, and more 
recent famous works [|TTj[ which essentially buried (at least, psychologically) non- 
supersymmetric unification. All analyses of the supersymmetric evolution of which 

1 am aware borrow the old technique used in the non-supersymmetric case |§. A 
new element which I would like to promote is the full exploitation of very specific 
features of the supersymmetric (3 functions UTTf . Among other virtues, this new ap- 
proach allows one to exactly solve (to all orders) the problem of the supersymmetric 
heavy thresholds. 

Supersymmetric grand unification, in the context of the realistic model-building, 
is briefly discussed in the second part. Here some simple numerical estimates are 
presented explaining the impact of various effects on a s [Mz)- 

2 Full Beta Function from One Loop of Pertur- 
bation Theory 

First, I would like to explain how the full multiloop (perturbative) f3 function can 
be obtained from a trivial one-loop calculation. The derivation presented below was 
reported over a decade ago |ff2|. Since this paper is hardly known to anybody it is 



instructive to reiterate the main points, the more so that we will need them later 
on, in discussing the threshold effects. 

Let us consider a generic non-Abelian gauge theory with the matter sector which 
includes (i) one chiral superfield in the adjoint representation of the gauge group (ii) 
additional matter superfields in arbitrary representations. For simplicity I will keep 
in mind that the gauge group G is SU(N) and assume that the additional matter 
is such that a mass term is possible for every superfield. Both assumptions can be 
lifted; moreover, the final expression we will obtain shortly is absolutely general. 

In more or less standard (and somewhat symbolic) notation the action has the 
form 

S = -X^Ti [ d 2 9d 4 xW 2 + J^Ti / d 2 9d 2 9d i x$e v <5> 
2qf, J 2<?g J 
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+ J J d 2 6d 2 6j2 S f ] e v S f + Qm° J d 2 9{S f f + h. c.) (1) 

where go is the bare coupling constant, $ is the superfield in the adjoint represen- 
tation, 

JV 2 -1 

$ = ^2 ® a T a , 

a=l 

Sf is the set of all other superfields, / is the "flavor" index, and m° is the bare mass 
term. Finally, T a are N x N traceless matrices, generators of SU(N). Let us note 
that the normalizations of the matter fields are different - $ includes the coupling 
constant while S/'s do not. As will become clear shortly, this is convenient. The 
bare Z factors of all fields Sf are set equal to unity. The lowest component of the 
superfield <3> will be denoted by ip. 

One may add all conceivable (super) Yukawa interactions between Sf] this does 
not change our argument. Perhaps, it is worth noting that in the absence of the ad- 
ditional Sf matter fields we actually deal with the extended N = 2 super symmetry. 

The model described by the action (p]) possesses JJ-flat directions - a system of 
degenerate classical vacua in the space of <3> fields , the so called valleys; the potential 
energy along the bottom of the valley vanishes. The valleys are parametrized by 
N — 1 chiral invariants, Tr $ 2 , Tr $ 3 , ... , Tr $ JV_1 . Alternatively, one can say that 
for any diagonal matrix <3> the D terms vanish. For instance, for 577(3) choosing 
if 3 = a complex constant and <p s = a complex constant and putting all other ip's to 
zero we get a point belonging to the bottom of the valley. 

There is only one point, namely the origin (ip a = for all a), where the full 
original gauge symmetry is unbroken. Any other point from the bottom of the 
valley corresponds to the spontaneous breaking of the symmetry. (This statement 
is valid, generally speaking, only in perturbation theory. Remember, however, that 
I am discussing only perturbation theory). Generically, SU(N) is broken down to 
{/(l)^ -1 . As a result, we have N — 1 massless "photons" with their superpartners, 
and N 2 — iV massive "W bosons", with their superpartners. These U W bosons" 
are analogs of the superheavy X and Y bosons of the standard theory of grand 
unification, sometimes called elephants. (Quite naturally - elephants live in the 
areas at the far side of the desert! All gauge bosons of the group G which acquire 
masses through the expectation value of the scalar field belonging to the adjoint 
representation of G will be referred to below as elephants.) Their masses need not 
be equal to each other. For example, for SU(3) 

Ml = |^ 3 | 2 , 

Mf )5 = i|V3y + ^| 2 , M 2 7 = i|V3y-^i 2 - (2) 

Now, we start our evolution from the normalization point /i equal to the ultra- 
violet cut off Mq, where the action coincides with the bare one, Eq. (|1|), and then 
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descend down from M to some low normalization point. M is assumed to be much 
larger than any other scale in the problem at hand. As we cross the mass thresholds 
we integrate out the corresponding massive fields, one after another. In this way we 
integrate out the elephants and all additional matter fields Sf. In the low-energy 
limit the only surviving fields are the N — 1 massless photons (and their superpart- 
ners), which are mutually neutral with respect to each other and do not interact. 
Since the residual gauge symmetry is U(l) N ~ l and there is no charged matter left 
below the last threshold, the gauge couplings do not run below the last threshold. 
There are N — 1 low-energy gauge couplings corresponding to N — 1 unbroken U(l) 
subgroups. These charges obtained after the completion of the evolution will be 
referred to as low-energy (frozen) charges a. They could be measured in a gedanken 
experiment, say, through the Coulomb interaction of heavy probe charged particles 
at large distances. 

We are free to choose any low-energy (frozen) charge and consider it as a func- 
tion of the ultraviolet cut off, the bare coupling constant ao = g\j^ and other 
parameters of the theory. Varying Mq and cto in a concerted way, ensuring that the 
low-energy quantities stay fixed, we find the M dependence of «o which is equiva- 
lent to the knowledge of the f3 function. For definiteness, one may concentrate on 
the charge corresponding to the "third photon" , A?, although this specific choice is 
of no importance for our purpose in this part. 

To connect a and «o one needs an explicit formula for the low-energy effective 
action, 

S = I ^ J d 2 6d 4 x W l W j + || J d 2 6d 2 6d 4 x ¥& . (3) 

This action includes only photons and the neutral adjoint matter. Correspondingly, 
the indices i,j run over the Cartan subalgebra (3 and 8 in the SU(3) example; in 
general they take rank(G) different values), and Z^ is the set of renormalization Z 
factors of the neutral adjoint matter fields evolved down "to the very end". Note 
that the kinetic term of the low-energy photons need not be diagonal in % and j. 

Limiting ourselves to the third photon gauge coupling (i = j — 3) one can 
trivially write the expression for 1/g 2 at one loop, 

Here T(Rf) is (one half of) the Dynkin index, 

Tr (T a T b ) = T(R)5 ab , 

T a are the group generators in the representation corresponding to the superfield 
Sf. For the adjoint representation of SU(N) the corresponding index T(adjoint) = 
T(G) = N. This fact is actually used in Eq. @; 2N in front of the first logarithm is 
actually 2T(G). For the fundamental representation of SU(N) the index T(fund) = 
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1/2 for each chiral superfield. It is worth reminding that one flavor requires two 
chiral superfields in this case, so that effectively T = 1 for each flavor. Finally, <p is 
a homogeneous function of the moduli whose concrete form depends on the group 
under consideration. In the SU(3) case for the third photon 



2/3 r i . •■> ; \ vn 1 (i 



This rather clumsy expression is nothing but a combination of masses ([2J). It is 
important to note that (p depends only on ip, not on (p^, and the only singular points 
(i.e. those where <p = 0) are 

ip 3 = or ip 3 = ±V3(p 8 . (5) 

At these points the pattern of the symmetry breaking is not generic; instead of 
SU(3) — > U{\) 2 we have SU(3) — > U(l) x SU{2), there are massless non-neutral 
gauge bosons, and the notion of the low-energy frozen action becomes ill-defined. 
We will stay away from the singular points. 

Now we come to a miracle. Although the low-energy coupling constant (|J) 
was obtained at one loop this expression is (perturbatively) exact to all orders - 
the gauge part of the effective action (Q) receives no correction at the level of two 
loops, three loops, and so on - at every finite order. This non-renormalization 
theorem stems from the fact that if the theory is fully regularized in the infrared 
domain there are no holomorphic anomalies Jl3|| , and the coefficient in front of W 2 



in the effective action must be an analytic function of g , m f and the moduli fields 



14] , |15l , |17ll . In the context of the string theory a similar observation was made in 
Ref. [TE|]. A more technical proof relying on specific properties of the background 
field technique for supergraphs was given in |15| (see also [|1(J) where the reader will 



also find a comprehensive list of earlier works in this direction. Note that the vacuum 
expectation values of the fields <p in all expressions above are those of the bare fields 
- the assertion of the holomorphic dependence refers to the bare parameters. Had 
the coefficient in front of W 2 been expressed in terms of the renormalized (low- 
energy) parameters, the holomorphicity would be lost. To be consistent, I should 
have marked the field dp by the subscript "0" ; I will do this occasionally, but not in 
every expression. The proliferation of indices makes them unreadable, so sometimes 
we will just keep this subscript in mind. 

It is very easy to see why higher corrections can not appear. Say, the two-loop 
term in the coefficient of W 2 , if it existed, should be proportional to l/(Re^ 2 ), 
since the imaginary part of g$ 2 (the 6 term) can not show up in the perturbative 
expansions. In which case there is no way to maintain the holomorphic dependence. 

Let us pause here and reiterate. The relation between the physical gauge coupling 
and the bare one presented above is (perturbatively) exact. The left-hand side is 
required to be independent of the ultraviolet cut off. This requirement determines 
the M dependence of g$ 2 and, hence, the full (3 function for the gauge coupling. 
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The (3 function is denned as 



(3 (a ) = da /din M . 

Naively, one might conclude that the one-loop expression for the physical charge (f|) 
yields the one-loop (3 function. Everybody knows that this conclusion is wrong, of 
course. So, there should be a subtlety. The question is "where ?". 

Within our formulation the change of g$ 2 under the variation of Mq should be 
adjusted in such a way as to keep all physical low-energy parameters fixed, not only 
g~ 2 . In particular, the physical (renormalized) masses of the matter fields and of 
the elephants must stay fixed; <p and m/ differ from <p Q and mj, however, by the 
corresponding Z factors, which, in turn, do depend on M . Therefore, differentiating 
the right-hand side of Eq. (|J) over In Mq it is necessary to differentiate ipo and m° 
as well. 

It is convenient to rewrite Eq. as follows 

1 "- Sir 2 if j 8ir 2 in j 

^ii^mi^). (6) 

vr 2 if y 8vr 2 m° f j y ' 

The sum of the first and the second brackets is holomorphic, each of them separately 
is not. The non-holomorphicity creeps in at this stage because the Z factors con- 
necting, say, rrif and m° are not holomorphic. The second line in Eq. (^|) reduces 
to 

1 (N . Za ^T{R f ) , \ . . 

— In h > y J ' InZt . (7) 

4tt \2tt a y 2n } J K ' 

In the absence of the additional matter fields supersymmetry extends to iV = 2, 
and the Z factor of the adjoint matter is Z = a^/a. In other words, (p is not 
renormalized in this case, and the j3 function is obviously one-loop. 

In deriving Eq. (0) I used the fact that the mass term of the additional matter, 




being the F term, is not renormalized ||19|| , and the entire renormalization of the 
mass of the /-th flavor comes from the corresponding Zf factor. With the additional 
matter fields included, the N = 2 supersymmetry degrades to iV = 1, and Eq. 
(0) generates the second and all higher-order terms in the (3 function. Indeed, 
differentiating ao over In Mq it is trivial to get 

1 2N ^T(Rf) 

JV ^T(R t ) \ 
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where 7 is the anomalous dimension of the adjoint matter, 7^ is that of Sf, 

7 = d In Z/d lnM 

and for simplicity it is assumed that there is no flavor mixing in Zy's due to (pos- 
sible) Yukawa interactions, which may be present, in principle, although I have not 
indicated them explicitly in Eq. (p]). Needless to say that this assumption is not 
crucial. The first and second lines in Eq. (|8]) come from the first and second lines 
in Eq. (El), respectively. Combining various pieces together we arrive at 



a 



0(a) = -7T 



NctY 
~2W 



f f 



(9) 



This is nothing but a particular case (occurring when one of the matter fields is in 
the adjoint representation) of the general NSVZ (3 function |TT| 



(3(a) = - 



a 
2^ 



T(G)a 
2vr 



3T(G)-^T(R i )(l- li ) 



(10) 



where in the expression above the sum on the right-hand side runs over all matter 
fields 



3 Where have the higher orders gone? 

Now, that I've explained how the multiloop f3 function appears, I am going to tell 
you that essentially it is never needed per se. Our task is calculating the gauge 
coupling; this calculation requires integration of the (3 function. But if we know the 
full expression for the gauge couplings beforehand why bother about getting the (3 
function and then integrating it back to find the couplings? 

If the evolution of the gauge couplings is complete and they are frozen, at the 
very "bottom" , one can represent the result as a pure one-loop expression provided 
the corresponding formulae are written in terms of the bare, not physical, threshold 
scales. Thus, if the bare values of the thresholds are given (as would be the case if 
a grand unified theory emerges, say, as a limit of string theory) we could just forget 
about the second and all higher loops altogether. 

Let me elucidate what I mean by this rather paradoxical statement. As an 
example, consider the SU(3) model described above, with one extra flavor in the 
fundamental representation (i.e. one chiral triplet superfield and one chiral an- 
titriplet). Since our task here is purely illustrative it is convenient to assume that 
the symmetry breaking 577(3) — > U(l) 2 takes place in two stages: first at a high 
scale SU(3) — > U(l) x 577(2) and then, at a somewhat lower scale, the remaining 
SU(2) is broken to U(l). In other words, instead of dealing with the general situa- 
tion (which can be addressed, of course) we will assume that |<£> 8 | ^> \(p 3 \- Among 
other simplifications this will allow us to neglect the off-diagonal term 

F(3) F (8) 

|Uf fill 
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in the low-energy effective action. Indeed, the coefficient in front of this term is 



3 / y 8 -(y 3 A/3) \ . 
32V3tt 2 V <P 8 + (^ 3 /V3) / ' 

in the limit at hand and is suppressed as |<^ 3 /<^ 8 |. 

Under the above pattern of symmetry breaking the elephants split into two 
groups - four very heavy elephants, with mass squared equal to 3|<^ 8 | 2 /4 emerging 
at the first stage and two not so heavy elephants with mass squared |v? 3 | 2 emerging 
at the second stage. There are four mass scales in the problem: the ultraviolet 
cut off Mq, the elephant masses and the (bare) matter mass m , with the following 
hierarchy 

m « |v9 3 | < \ip 8 \ < M . 

As we descend from Mq, we pass the first elephant threshold where the gauge group 
SU(3) is broken down to U{\) x SU{2). The matter triplet becomes an 577(2) 
doublet plus a singlet coupled only to the eighth photon. Then, below the second 
elephant threshold, only the Abelian subgroup survives; we have two photons and 
three matter fields. The charges of these three matter fields with respect to these 
two photons are 

/l 1 \ 

I - , — - , 1 for the third photon, 

and 



(— — 

\2 v / 3 ' 2v^3 ' V3 



— -= ] for the eighth photon. (11) 



The gauge couplings a-$ and as continue to evolve due to the contribution of the 
matter fields. 

When we cross, in our descent, the first elephant threshold we get two separate 
(3 functions governing the evolution of as and the SU(2) gauge coupling. Then, at 
the second elephant threshold the SU{2) gauge coupling becomes a%. Below the 
matter threshold the evolution of the gauge coupling constants stops, and we arrive 
to what we call frozen couplings. 

The final answer for the low-energy (frozen) coupling constants is known exactly, 
as was demonstrated above. The all-order result for a^ was given in Eq. (|]). For 
the sake of convenience we reproduce it here again, along with the answer for a§, 

1 _ 1 6_ ln M i 1 in M ° 

«s « 2vr (^3^/2) 2f m° ; K ' 
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1/a 



l/a 3 



l/a 8 




Zn(/x) 








Figure 1: Evolution of the gauge couplings: from unification to freezing. The straight 
lines represent the one-loop evolution with the fake (bare) thresholds. The actual 
evolution is given by the smooth curves. 

In the difference of the constants nearly everything drops out, 



I emphasize that it is the ratio of the fields' bare expectation values that enters 
here. This result looks remarkably simple! One could superficially interpret it is 
follows. At the first (heaviest) elephant threshold, \/3ipl/2, the couplings 1/a are 
unified. They diverge immediately, running linearly (in the log scale), according to 
the one-loop formula, all the way down. The slope of the linear running changes 
abruptly at the elephant thresholds. From 5 above the heaviest elephant threshold 
it changes to 3 for the SU(2) coupling and —1 for the eighth photon. Below the 
second elephant threshold, <^g, where we deal with the U(l) couplings, they run in 
parallel, in accordance with the charge assignment ([O]), until both couplings reach 
one and the same threshold at mo, where they freeze at once (Fig. 1). This simple 
interpretation is operationally correct, but this is deceptive correctness, of course. 

After briefly reflecting one gets a second thought. First of all, it is absolutely 
certain that the (3 functions are multiloop. Therefore, the running law is non-linear, 
and this non-linearity is different for l/a^ and l/a%, in particular, in the interval 
between the two elephant thresholds since the Z factors of the SU(2) doublet and 




(13) 
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singlet in this domain are definitely different. Second, for the same reason the matter 
threshold can not be at one and the same place for the third and the eighth photons 
- although we start from one and the same parameter m the physical masses of 
the three fermions we get after SU(3) is broken down to U(l) 2 are different. There- 
fore, the constants must freeze out at different scales. Finally, v^Vo/^ and are 
not physical threshold scales as well; these are bare parameters, the corresponding 
physical thresholds differ by the Z factors. 

And yet, all these complications (the difference between the physical and bare 
thresholds and the higher orders in the f3 functions) conspire together and kill each 
other producing an effective picture of the one-loop evolution with the slopes jump- 
ing at fictitious thresholds. The actual evolution can be obtained from the fictitious 
straight lines of Fig. 1 by slightly distorting the scale in the horizontal direction. 
The distortion is logarithmic, i.e. rather weak (provided all couplings enter in the 
weak regime everywhere below Mo) and is scale (and pattern) dependent. In the 
SU(3) and SU(2) parts we pull the thresholds to the left, while in the U(l) parts 
we pull them to the right, by different amounts. The distortions of the /i scale are 
determined by the Z factors which, although different from unity, are not strongly 
different. Under realistic conditions they may be 1.5 or 2 or so. 

Thus, we witness another little miracle. Although the derivation presented above 
refers to a particular case, actually the message is general - (i) use the bare values of 
the threshold parameters abstracted from the bare Lagrangian God-given at some 
high (say, Planck) scale; (ii) forget about the multiloop effects, just run and match 
the gauge couplings according to the most naive one-loop formula (the so-called d 
function approximation). Then, at the very end, when the evolution is finished, you 
will get the correct full results for the low-energy gauge constants. 

In particular, let us draw attention to the fact that the additional matter pro- 
duces no impact whatsoever in the difference of the couplings, Eq. (|13|) . If we 
worked in terms of the physical thresholds, we would start from different physical 
masses for the three species of the additional matter. The Z factors will also be 
different. But the product of m's times Z's will be the same. That's how the ad- 
ditional matter contribution will disappear in the difference of the gauge couplings, 
although it will be much harder to see this disappearance keeping in mind that the 
Z factors are usually calculated approximately, not exactly. 

This feature is general. Say, in a realistic SU(5) model all SU(5) multiplets 
provided by a common mass parameter are irrelevant for the calculation of the 
coupling differences. 

What will happen if we introduce physical threshold scales - the "GUT" scales 
V3(p 8 /2, instead of the bare scales? Then 



plus logarithm of the Z factor describing the evolution of the field in the interval 
between the first and the second elephant threshold. If we transfer this logarithm 
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from the right-hand side to the left-hand side of Eq. (p~3|) , on the left-hand side we 
will get a specific combination which is actually nothing else than the Wilsonean 
couplings |T5|]. The issue of the thresholds will be discussed in more detail later on. 



4 What if the evolution is not completed? 

The example discussed above is very elegant but, unfortunately, is not very practical 
since usually the symmetry breaking pattern is such that not all gauge bosons acquire 
masses, so it is impossible to descend down, below the " freezing point". If the 
freezing points existed, the scale where they (the freezing points) occur is established 
as a result of a natural flow, and is different in different subgroups. This variation of 
scales neutralizes the effect of the second and higher-order terms in the (3 function. 

In the case there is no freezing point, or we measure the running couplings above 
the freezing points, we usually measure them at a scale /i which is one and the same 
for all subgroups. Then the impact of the second and higher-order terms in the j3 
function is nonvanishing. To get a better idea of what is going on we will consider 
another pedagogical example, the SU(3) model we dwell upon above, but this time 
it will be assumed that the symmetry breaking pattern is SU(3) — > SU(2) x U(l). 
Moreover, for the time being we will discard the additional matter, and add the 
superpotential term to the $ superfield 

W = A (™o$ a $ a + Ao<U c $ a $ 6 $ c ) , (14) 



where d a b c are the d symbols of SU(3). The superpotential ([14]) destroys the valley 



and fixes the vacuum value of the $ field. One of the possible solutions is 

8 2m 

this solution is characterized by the property of vanishing if 3 which ensures that 
the SU(2) subgroup remains unbroken. The bare mass squared of the fourth, fifth, 
sixth and seventh gauge bosons, the elephants, is 



(Mf 7 ) 2 = (3|^| 2 )/4 



they eat up the corresponding $ superfields. The first three gauge bosons remain 
massless. There are two thresholds - one associated with the elephants, another one 
with the fields $ 1,2 ' 3 . The latter are massive because of the superpotential (|i~4l) and, 
in particular, the mass term of the $ superfield. The bare mass of $ 1,2 ' 3 is 



K> 3 ) : 



1 2 



Let us assume that the elephant threshold is heavier than that of $ 1,2 ' 3 , which 
is natural if the coupling constant A is weak. Then below the lowest threshold 
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associated with (J) 1 ' 2 ' 3 we are left with SU(2) supersymmetric gluodynamics, with 
no matter fields, plus the eighth photon and no charged matter. So, the eighth 
gauge coupling is frozen; the SU{2) coupling, however, evolves till the very end - 
this subgroup, being unbroken, has no natural freezing point. 

Thus, far above the thresholds we start from the unified gauge coupling chq- 
Far below the two thresholds we have the evolving coupling, asu(2) ; and the frozen 
coupling, otu(i). Since the SU(2) coupling evolves we have to distinguish between 
the Wilsonean constant - it will be marked by braces - and that appearing in the 



c-number functional T [15]. The Wilsonean constant {!/«} is renormalized only at 



one loop, while 1/a is renormalized to all orders. The exact relation between 
and 1/a is explicitly known |TJ 



Let us choose the normalization point // far below the lowest threshold, [i <C 
M]_3. We start from the following one- loop formula for the Wilsonean couplings 
(which is perturbatively exact) 

f_ar I 2,_ 61n Mo_ 31n Mo + 
l«su(2)Wj «o A* (V3^g/2) 

, M , M 
21n— + ln^ , 15 

m (V3^/2)' 1 ; 

and 

2vr 2tt , M 
= 61n^ . 16 

«u(i) "o (V3<^o/2) 

The first line in Eq. (|T5j ) reflects the contribution of the gauge bosons, and the fact 
that the fourth to seventh become elephants and freeze out. The second line in Eq. 
(|I5|) is due to the <3? fields. The first term in the second line comes from $ 1,2 ' 3 , the 
second from those $'s ($ 4 to $ 7 ) which are eaten up by the elephants. We know 
already that Eqs. flTJD and ([TBI) are valid to all orders. 

The expressions presented above are not so easy to grasp; they may even generate 
a couple of perplexing questions. To get a transparent interpretation of this result 
two further steps are in order. First, we will pass from the Wilsonean coupling 
|i/q!su(2)(a*)| to the regular one, l/asu(2)(A i )- Second, it is convenient to consider 
the difference between l/asu(2)(/-0 an d l/au(i)> keeping in mind that the latter is 
frozen below the elephant threshold. Therefore, «u(i) i n the problem at hand is 
nothing else than ckgut- 

In pure SU(2) supersymmetric gluodynamics the first transition is realized as 
follows [0 

f^l = ^ + 21n^, (17) 
L«SU(2)(j"J J asu(2)(/f) "0 

where the factor 2 in front of the logarithm is actually T(G) for G = SU(2). In 
addition, we have to introduce the Z factors describing the evolution of the $ matter. 
Let Z(M — > Mqut) express the evolution from M down to the physical GUT scale 
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M 4 _7 = Mqut, while Z(M — > Mi_ 3 ) expresses the evolution from M down to the 
physical m threshold (for to $ 3 ). Then, with our definitions, m = Z(M — > 

Mi_ 3 )m and v^o/ 2 = M GVT [Z{M -> M G uT)]~ 1/2 a2 /2 aGUT- Subtracting Eqs. 
(P^D and ( |16D from each other and using the physical thresholds we readily get 

fr = Ji!— 61n M ^ , ; , + 21n - M °" T „ (18) 

«SU(2)(/i) «GUT /^_ogu2_V mZ(M GU T^Mi_ 



\ y a SU(2)(M) 

where the factor 

Z(M GVT -> Mx_ 3 ) = Z(M -> M 1 _ 3 )/^(M -> M GUT ) 

expresses the evolution from the elephant (GUT) threshold to the m threshold; thus, 
the product mZ (Mqut — ► Mis) is nothing else than the mass parameter Mi_ 3 at 
the GUT scale (rather than the physical mass of $ 1 ~ 3 ). 

Equation ( |18D has various appealing features. In particular, the parameters M 
and a referring to the initial scale, lying far above the GUT scale, drop out, as they 
should. 

Now, the last exercise aimed at making our toy model of grand unification a little 
bit more realistic. Let us add an extra matter superfield, say, in the fundamental 
representation of SU(3) (actually we will need two chiral superfields - one triplet 
and another antitriplet), with no mass term. What happens then? 

A qualitative impact is quite obvious - the U(l) coupling constant now runs 
below the GUT threshold, and never freezes. Therefore, it is necessary to intro- 
duce now the Wilsonean coupling for U(l) as well. Quantitatively, an extra term 
appearing on the right-hand side of Eqs. ([15]) and (|T6| ) is 

In , 

one and the same in both expressions. Simultaneously, 27r/a!u(i) m Eq. ( |16|) should 
be substituted by {27r/au(i)}. Transition from the Wilsonean couplings to the reg- 
ular ones is now to be carried out in both subgroups, SU(2) and U(l). In the first 
case, instead of Eq. (O) we have 



= — ~~~r^ + 2 ln + ln Z ( M o - n d.m.) , (19) 

«su(2)CwJ asu(2)(yw) "o 



1 



where d.m. stands for the doublet matter (s.m. below stands for singlet matter, 
i.e. the matter field from the original triplet which does not belong to the SU(2) 
doublet). 

Furthermore, for the U(l) subgroup 

1 1 +- In Z(M -> fx;d.m.) + - In Z(M ^ fi; s.m.). (20) 



1 



«u(i)(^) J q u(i)(aO 3 ' " 3 
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From M down to Mgut the Z factor is common, below the GUT scale it splits. The 
SU(2) doublet interacts with both, the SU(2) gauge bosons, and the eighth photon 
(see Eq. (0)). The SU(2) singlet interacts only with the photon. 

Assembling all pieces together we observe that after the additional matter is 
included, the following relation emerges 

2tt 2tt 



«SU(2)(/i) «u(i)(aO 
-6 In 775- + 2 In — t—— — — - + m ■ 



"GUT 
a SU(2) O) 



V3 mZ (Mgut -> Mi_ 3 ) /^(M GUT //; d.m.) 



11 /x[Z(M GUT - //; rf.m.)] V3[ Z (M GUT - /i; s .ro.)] 2/3 ' 1 j 

where each term in each line has a very clear physical interpretation. The three terms 
in the second line are due to the SU(2) gauge bosons, the heavy mass threshold of 
the $ fields; and the additional matter; the third line is the evolution of the U(l) 
gauge coupling. 

Our toy example is now generic enough in the sense that it includes all elements 
one may encounter in the realistic models. Therefore, at this stage we are able to 
formulate general conclusions. We start from the discussion of the mass thresholds. 



5 Mass thresholds 

The proper treatment of the heavy threshold effects in the evolution of the gauge 
couplings is crucial in establishing relations between the low-energy couplings fol- 
lowing from the fact of unification at a higher scale. At the one-loop level there is 
no problem - the so called i? function (or run-and-match) approximation is appli- 
cable. According to this approximation, as /i evolves from Mgut downwards, one 
just abruptly changes the first coefficient of the (3 function (i.e. the coefficient in 
front of the logarithm) at the physical value of the corresponding mass. Say, at 
the $ threshold, above M$ one includes in the coefficient the $ contribution while 
immediately below M$ one excludes the $ fields altogether. At the one-loop level 
the $ function prescription is certainly correct. 

The question arises beyond the leading log approximation where the above 
prescription is certainly incorrect. Following a tradition established in the pre- 
supersymmetric era the same i? function approximation is usually applied at two 
loops |2lJ. Even though numerically this may not be such a bad idea, it is interest- 



ing to ask what is the correct way of taking the heavy mass thresholds in the gauge 
couplings into account. 

An idea which immediately comes to one's mind is smoothing the d functions 
in some way [ 22| . The smoothing procedure is rather ambiguous, it depends on 
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conventions and, generally speaking, has no direct physical meaning. This happens 
because the gauge couplings by themselves are not physical observables. The defi- 
nition we follow - through the coefficient in front of W 2 in the effective action - is 
directly related to observables (with the power accuracy) only provided we are far 
above the threshold or far below it. Fortunately, with respect to the heavy thresh- 
olds, we are always far below. If we calculate at low energies the effect due to the 
heavy fields (which are integrated out since they can appear only in virtual loops) 
in the full theory it is all we need. And, as was shown above, that's exactly what 
we do. 

Let us examine Eq. ( |21| ) (the same conclusions follow from the analysis of Sect. 
3). The heavy threshold effect in this expression is represented by the second term 
in the second line, 

2 j n ^GUT 

11 mZ(M GUT -> Mi_ 3 ) ' 
If the physical mass of the fields $ 1-3 coincides with Mqut (i-e. with the elephant 
masses), then the argument of the logarithm is unity - there is no heavy threshold 
correction at all. This is quite natural and uninteresting. What is more interesting 
is the case when the physical mass of the fields $ 1-3 is less than Mqut- Then the 
threshold correction is non- vanishing. In the $ function approximation the threshold 
correction will reduce to 

r, , ^GUT 

2 In 

m 

where m is the physical mass of the fields $ 1-3 (the same as Mis). 

Now, the correct answer, taking into account all orders, replaces the physical 
mass m by itiZ^Mgvt — * M1-3), or the corresponding mass parameter at the GUT 
scale. If our research starts from model-building and goes in the direction of phe- 
nomenology, the short distance (GUT) value of the mass parameter is primary, and 
no further efforts are required. Case closed. If, however, we start from phenomenol- 
ogy, and operate with the physical masses of the particles, then to properly include 
the heavy threshold effects beyond one loop (i.e. beyond the •§ function approx- 
imation) one has to calculate the corresponding Z factors. Unfortunately, the Z 
factors are not calculable to all orders. Note, however, that in the two-loop analysis 
of the gauge couplings the Z factors appear only at one-loop (i.e. in the leading log 
approximation). In other words, we reduce the level of complexity by one order. In 
particular, if there are heavy threshold corrections in the Z factors themselves ( they 
occur more rarely than in the gauge couplings, if at all) it is legitimate to account 
for them in the d function approximation. This will provide us with the two-loop 
evolution of the gauge couplings properly accounting for the heavy thresholds. 

6 General lessons and master formula 

In order to analyze the evolution of the gauge couplings from M G ut down to some 
low scale, first one starts from the standard one loop expressions written in the most 
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primitive "run and match" approximation, using the physical values of the mass 
parameters in the logarithms. By the mass parameters I mean the threshold masses, 
^gut> and the normalization point /i. Then, in the logarithms corresponding to 
the matter field contribution, substitute the mass parameter by the mass parameter 
multiplied by the corresponding Z factor. In the logarithms corresponding to the 
gauge field contribution, substitute the mass parameter by the mass parameter 
multiplied by the corresponding Z 1 ^ 3 factor. For the gauge fields the Z factor is 
l/a. 

The above procedure fully specifies the evolution process. It is perturbatively 
exact to all orders, and takes into account the heavy mass thresholds exactly pro- 
vided that supersymmetry holds, i.e. we stay above the scale of the supersymmetry 
breaking. 

(In practical calculations one usually limits oneself to two-loop accuracy. Even 
for that limited purpose our master formula is extremely useful. First, it spares one 
the necessity of tabulating the two-loop coefficients of the (3 function, a rather cum- 
bersome task by itself. Only one loop coefficients enter the game, and all one loop 
coefficients are simple numbers having a geometric, and group-theoretic, meaning. 
Second, it provides with a consistent and exhaustive treatment of the mass thresh- 
olds, as was mentioned above. Third, numerical integration of the renormalization 
group equations becomes unnecessary. The final result is expressed in terms of the 
Z factors. In the two-loop analysis the latter enter only at one loop level, or, more 
exactly, in the leading log approximation, which is equivalent to one loop.) 

In the text book case of the SU(5) grand unification |23[ the master formula is 



2vr 2tt „, Mqut 



, - 9 In — + E 

«GUT ( "gut A 1 gen L 



in 1 — In 1 — in 



liZ qL 2 fiZ UR 2 fiZ 



DR 



+ 



Mgut , M GVT 

+ 31ll liir + ln ]Sfr' (22) 
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GUT 
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M GUT 4 M GUT 

2 111 7^— + - 111 frr- 



(24) 



Here I made various simplifying assumptions (which, in principle, can be easily 
lifted). First, it is assumed, of course, that \i is larger than the superpartners' 
masses - our master formula is valid only as long as supersymmetry holds. The 

(2) 

masses of the two Higgs doublets of the model at hand, Mjj are also assumed to 

(2) 

lie above the supersymmetry breaking scale, and above /x. Let us stress that M H 
in the expressions above is not the physical mass, but, rather the corresponding 
mass term at the GUT scale. At the same time the SU(2) breaking scale is taken to 



be below fi - the SU(3) x 5*7(2) x U(l) symmetry is manifest in Eqs. (||) - (||). 
Moreover, the Yukawa interactions are neglected so that all three matter generations 
have identical Z factors. This means that in the approximation accepted there are 
only five different matter Z factors - Z q i corresponding to the left-handed quark 
doublet, Z UR corresponding to the right-handed up-quark, Z DR corresponding to the 
right-handed down-quark, Zu, corresponding to the left-handed lepton doublet, and, 
finally, Z eR corresponding to the right-handed electron (muon, r). Mjj stands for 
the mass term of the triplet superheavy Higgs particles (from 5 and 5), as it is seen 
at the GUT scale. Finally, M$ is the mass term of the chiral superfields from the 
24-plet at the GUT scale. The last line in each of the expressions above represents 
the GUT scale threshold terms, plus the doublet Higgs contribution in Eqs. (0), 

All the simplifying assumptions listed above can be lifted in an obvious way. 



7 Z factors 

The master formula above expresses the gauge couplings, to all orders, in terms 
of the Z factors. Unfortunately, the latter are not calculable to all orders. The 
Z factors appear as the coefficients in front of the D terms, and the theorems of 
the holomorphic dependence on the moduli and the (inverse) gauge couplings are 



not valid p4| . Without the power of holomorphicity one can not establish the 
all-order results for the Z factors and has to resort to the old-fashioned order-by- 
order perturbative calculations. The anomalous dimensions of the matter fields 
were calculated up to three loops |25|. The higher-order terms in 7's are definition- 



dependent, however, starting from the second loop. The work of bringing the results 



in line with the definition where the NSVZ j3 function is valid is under way now [26| . 
After this work is done one will be able to calculate the gauge coupling evolution 
up to fourth order. 

As was mentioned more than once, for the two-loop analysis of the gauge cou- 
plings, it is sufficient to have the Z factors in the leading log approximation, which 
is unambiguous and scheme-independent. Actually, the one loop expressions can be 
obtained with no calculations, from general arguments alone. 
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Figure 2: One- loop graph describing the matter field renormalization due to the 
exchange of the gauge field. 

The relevant one- loop diagram is presented on Fig. 2. Let us assume at first 
that we have only one chiral matter superfield, in the adjoint representation, with 
no self-interaction. Then the theory has N = 2 supersymmetry, and the (3 function 
is well-known to be one-loop. From Eq. (|T0| ) we then immediately read off that 



7 = —N— . 

From the structure of the graph depicted on Fig. 2 it is clear that the factor N 
above is nothing else than (^(adj) where C 2 is the quadratic Casimir operator, 

(T a T a ) rcpR = C 2 (R)l, 

related to the Dynkin index in a standard way, 

dim (adj) 



C 2 (R) = T(R)- 



dim (R) 



At the one loop order the diagram of Fig. 2 is the only relevant graph. It is obvious 
then that in the general case of the matter field in the arbitrary representation 
of the gauge group, C 2 (adj) is merely substituted by C 2 (R), irrespective of how 
many matter fields we have, whether or not there are Yukawa interactions (self- 
interactions), and so on. We conclude that the one- loop gauge contribution to Z 
is 

Z = l-C 2 (R)^\n^. (25) 

It is clear that the gauge-interaction contribution is negative, tending to make Z J s 
less than unity. 

The renormalization group improvement of expression fl2"5]) (summation of the 
leading logarithms) is trivial, 

Mi) • (26) 

where 

b = 3T{G)-Y,T{Ri) 
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is the first coefficient in the (3 function. 

Similar "back-of-the-envelope" arguments yield the (super)Yukawa vertices' con- 
tribution to Z's, which, on the contrary, tends to make the Z factors larger than 
unity Indeed, let us start again from the case of adjoint matter, this time assuming 
that there are three matter fields, an d If the Yukawa coupling 



AW = ^f abc ^ b 2 % 



(27) 



is added to the superpotential the theory becomes N = 4 supersymmetric, the gauge 
coupling does not run, and the anomalous dimensions of all $ fields vanish. This 
means that the renormalization of, say, the $1 field due to the gauge boson (Fig. 
2) is canceled by a similar diagram with the gauge boson line replaced by that of 
the <3> fields, and the gauge vertices replaced by the (super)Yukawa vertices. We 
immediately conclude then that 



J Yukawa 



1 + A^ln^ 

7T fl 



(28) 



Getting the Yukawa contributions from Eq. yzq) for other matter fields and different 
structure of the Yukawa vertices is a mere substitution of obvious color factors. 

It is convenient to collect here the leading-logarithm expressions for the five Z 
factors entering our master formula in the SU(5) case, 



ZsU(3) X Z su(2 ) X Z 1 



1/9 



^UR ~~ 

Zdr 



Zsv(3) x Z 1 



16/9 



J SU(3) 



x Z\ 



1/9 



and, finally, 



where 



^SU(3) 



— ^SU(2) X Z\ 

Z P - ^ 



>eR 



z\ 



-1 8/9 



J SU(2) 



a 2 (fi) 



-3/2 



&GVT 

a x (fi) 



-1/22 



(29) 



(30) 



and the Yukawa interactions are neglected, for simplicity. The Yukawa interactions 
might be important for the t quark, but even in this case, numerically, their effect 
is very modest. Moreover, it is assumed that the heavy thresholds coincide with the 
GUT scale, while the doublet Higgs particles are lighter than fi. It is quite clear 
how to amend this expression if /1 is smaller than the doublet Higgs masses. 
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8 Numerical analysis: bird's eye view 



Now, we will discuss the numerical situation, in gross features, leaving aside fine 
details. Our task is merely to get a general idea. To this end I will neglect the low- 
energy thresholds due to supersymmetry breaking, assuming that all superpartners 
have masses in the ballpark of 100 to 200 GeV. Making some of them significantly 
heavier (say, ~ 1 TeV) may change our conclusions. This is a different topic, how- 
ever, which I will not touch upon, referring the reader to the very rich literature 
devoted to the issue of low-energy thresholds. Two other possibilities - new physics 
in the middle of the great desert, and Planck physics corrections at the GUT scale - 
will not be considered as well. The first topic was covered in the talks of B. Brach- 
machari, E. Keith and some others, while the second topic was extensively reviewed 
by Prof. P. Nath. 

To begin with, let us truncate the expressions (p^ ) - (|24D at one loop and dis- 
regard the heavy thresholds (i.e. assume that the masses of the Higgses from the 
24-plet and triplets from the quintets coincide with Mqut)- Then for the low-energy 
value of a 3 we get 

2?r _ 12 2tt 5 2tt 

— « i (31J 

as 7 o>2 7 Oi\ 

where by low energy I mean that the corresponding normalization point is somewhere 
in the ballpark of 100 to 200 GeV, see above. The value of 2n/a2 is close to 186 
and that of 2ir/ai is close to 371. Where exactly the normalization point lies affects 
the numbers above at the level ~ one unit; a similar uncertainty is associated with 
the low-energy thresholds (provided that the sparticle masses are ~ 100 GeV). The 
effects I am hunting for now constitute seven units - this is the impact of the 
difference between a s (M z ) = 0.11 and a s (M z ) = 0.125. 

Using Eq. ([H]) we conclude that 27r/a 3 « 54 (and 27t/o;gut ~ 153). If the 
second loop shifted this number up by three units, to 57, I would be very happy. The 
above value then would nicely match what is expected from the low-energy QCD 
phenomenology, a s [Mz) ~ 0.11 |J. Unfortunately, the higher-order corrections 
work just in the opposite direction shifting 277/0:3 down by three units. 

By inspecting Eqs. ([22] ) - (|24| ) and ( |29| ) we immediately observe that, by far, 
the most significant numerical effect of higher loops comes from the factors 

\ 1/3 / \ 1/3 

<*M) VM/i)/ 

in the logarithms in Eqs. (|22|), (f2~3|) and from Zsu(3)- The factors Zi and Zsu(2) 
in the standard version of the SU(5) grand unification are too close to unity to be 
of importance in my back-of-the-envelope analysis - they are buried in the noise of 
other insignificant corrections. Keeping only the important factors we get, instead 
ofEq. fl3J, 

2tt _ 12 2tt 5 2?r 

«3 7 «2 7 OJi 
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In the second line we can use the values of the constants in the one-loop approxima- 
tion; moreover, Zsv(3) is given by Eq. (|30|) . The most significant contribution, —3, 
comes from the first term in the second line. The second and the third terms are 
close to 0.6 and nearly compensate each other. As a result, the second loop shifts the 
prediction for 277/0:3 away from the right value, from 54 to 51 (i.e. a s (Mz) ~ 0.123). 

Since we are so remarkably close to the desired number it is quite natural to 
think that the residual discrepancy is due to some corrections that are unaccounted 
for so far. The idea which immediately comes to one's mind is including the Yukawa 
interactions in order to change the values of the Z factors. We - Ian Kogan and 
myself - played with this idea for some time. To achieve success in this way one 
has to make an appropriate Z factor of order 100 or larger (InZ ~ 5). It is quite 
obvious that in the minimal model even the strongest Yukawa interaction, that of 
the t quark, is far too weak to ensure such a large value of Z. Only if some Yukawa 
constant approaches the Landau pole can we expect to get Z that large. 

Another way out may be associated with the heavy particles' thresholds. Gen- 
erally speaking, the masses of the Higgs particles from the 24-plet and the triplets 
from the quintets need not exactly coincide with Mqut, see the second lines in Eqs. 
(p2|) - fl2~4]). Again, it is clear that to shift the prediction for 2ir/a by several units 
the logarithm of the mass ratio has to be several units itself. In other words, the 
superheavy Higgs particles have to be hundred times lighter than Mqut- In this 
case there is the menace of a fast proton decay due to dimension five operators [27j. 



In more details the scenario was investigated in Refs. [28, 29, 30]. It was shown 



that in the SU(5) grand unification this solution does not work, unfortunately. The 



proton decays too fast [28, ^9|. In the SO(10) grand unification, however, the super- 
heavy Higgs thresholds do bring the value of a s (Mz) down to 0.11 without making 
the proton lifetime unacceptably short p0| . Moreover, the prediction for the proton 
lifetime returns to the range accessible experimentally, which makes the whole story 
exciting. 



9 Conclusions 

• If the primary "high energy" gauge group G is spontaneously broken, G — > 
G\ X Gi x then the low-energy couplings l/oi, I/0J2, ••• are given by a master 
formula to all orders, including mass thresholds. 

• If supersymmetric mass parameters are known at the GUT scale (rather than 
the physical mass parameters) then full inclusion of the threshold effects trivializes. 

• If the evolution is complete (i.e. all GVs are f/(l)'s, and the low-energy cou- 
plings are frozen, then the all-order result for I/014 looks as if it was one-loop, with 
faked (bare) thresholds. 
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• The center of gravity of the multiloop analyses is shifted from the gauge (3 
functions to the Z factors of the matter fields. The complexity (number of loops) is 
reduced by one level. In the minimal SU(5) grand unification the dominant two-loop 
correction comes from the one-loop Z factor of the gluon (gluino) field. It is of the 
right magnitude, but, unfortunately, its sign is "unfavorable", so that the value of 
a s (M z ) becomes unacceptably high. 

10 Acknowledgments 

I would like to thank I. Kogan, L. Roszkowski, and A. Vainshtein for stimulating 
discussions. Assistance of B. Chibisov with plots is gratefully acknowledged. This 
work was supported in part by DOE under the grant number DE-FG02-94ER40823. 

References 

[1] The paper of Gol'fand and Likhtman was published in April 1971 (JETP Lett. 
13, 323). The work of Wess and Zumino, who came to supersymmetry from 
a totally different direction and were unaware of Gol'fand and Likhtman, ap- 
peared in the end of 1973 (Nucl. Phys. B70 (1974) 39). 

[2] M. Shifman, Mod. Phys. Lett. A10 (1995) 605. 

[3] M. Shifman, Int. Journ. Mod. Phys. All (1996) 3195. 

[4] A. Kataev, A. Kotikov, G. Parente and A. Sidorov, Next-to-next-to-leading 
Order QCD Analysis of the CCFR Data for xF% and F2 Structure Functions 
of the Deep Inelastic Neutrino Nucleon Scattering, Preprint INR-0919/96 |[hcp- 
| ph/ 9605361 . 

[5] The most uncompromising point of view was seemingly formulated by G. Kane, 
who summarized existing evidence in favor of a supersymmetric model with 
the relaxed constraint of the universal gaugino (supersymmetry breaking) mass 
terms. Incidentally, the non-universality of the gaugino mass terms also helps 
get the right value of a s (Mz) in the most naive version of grand unification, 
see e.g. L. Roszkowski and M. Shifman, Phys. Rev. D53 (1996) 404. 

[6] H. Georgi and S.L. Glashow, Phys. Rev. Lett. 32 (1974) 438. 

[7] See e.g. the recent review papers 

P. Langacker, Tests of the Standard Model and Searches for New Physics, in 
Precision Tests of the Standard Electroweak Model, Ed. P. Langacker (World 
Scientific, Singapore, 1994); 

G. Kane, C. Kolda, L. Roszkowski, and J. Wells, Phys. Rev. D49 (1994) 6173; 
M. Carena, Unification of Couplings in the MSSM, Proc. Int. Europhysics Conf. 



22 



on High Energy Physics, Brussels, 1995, Eds. J. Lemonne, C. Vander Velde and 
F. Verbeure (World Scientific, Singapore, 1996), page 481 |[hep-ph/9603370|| ; 
and references therein. 



[8] H. Georgi, H. Quinn, and S. Weinberg, Phys. Rev. Lett. 33 (1974) 451. 

[9] M. Einhorn and D.R.T. Jones, Nucl. Phys. B196 (1982) 475; 
W. Marciano and G. Senjanovic, Phys. Rev. D25 (1982) 3092. 

[10] J. Ellis, S. Kelley and D. Nanopoulos, Phys. Lett. B260 (1991) 131; 
P. Langacker and M.X. Luo, Phys. Rev. D44 (1991) 817; 
U. Amaldi, W. de Boer and H. Furstenau, Phys. Lett. B260 (1991) 447; 
F. Anselmo, L. Cifarelli, A. Peterman and A. Zichichi, Nuov. Cim. A104 (1991) 
1817. 

[11] V. Novikov, M. Shifman, A. Vainshtein, and V. Zakharov, Nucl. Phys. B229 
(1983) 381; Phys. Lett. B166 (1986) 329. 

[12] A. Vainshtein, V. Zakharov and M. Shifman, Yad. Fiz. 43 (1986) 1596 [Sov. 
J. Nucl. Phys. 43 (1986) 1028], Sect. 3; see also Pis'ma ZhETF 42 (1985) 182 
[JETP Lett. 42 (1985) 224]. 

[13] M. Shifman and A. Vainshtein, Nucl. Phys. B359 (1991) 571. 

[14] E. Cremmer, S. Ferrara, L. Girardello and A. van Proeyen, Nucl. Phys. B212 
(1983) 413. 

[15] M. Shifman and A. Vainshtein, Nucl. Phys. B277 (1986) 456. 
[16] M. Shifman and A. Vainshtein, Nucl. Phys. B296 (1988) 445. 
[17] N. Seiberg, Phys. Lett. B318 (1993) 469. 

[18] L. Dixon, V. Kaplunovsky and J. Louis, Nucl. Phys. B355 (1991) 649. 

[19] M. Grisaru, M. Rocek and W. Siegel, Nucl. Phys. B159 (1979 429. 

[20] If all mater fields we have in the model belong to the fundamental representation 
of the gauge group it is possible to prove that the NSVZ (3 function is exact 
also including nonperturbative effects. The proof can be abstracted from Ref. 
|T6| . In this work the gaugino condensate (AA) was exactly calculated. Since 
(AA) is in the same superfield as the anomaly of the trace of the stress tensor 
(AA) has zero anomalous dimension; if it is known exactly the full /3 function 
immediately follows. 



23 



[21] The literature on this subject is so vast I can not even scratch the top of the 
iceberg. Some recent publications are: 

V. Barger, M.S. Berger, and P. Ohnmann, Phys. Rev. D47 (1993) 1093; 
D. Castano, E.J. Piard, and P. Ramond, Phys. Rev. D49 (1994) 4882; 

G. L. Kane, C. Kolda, L. Roszkowski, and J. Wells, Phys. Rev. D49 (1994) 
6173; 

M. Carena et al, Phys. Lett. B317 (1993) 346; 

L. Clavelli and PW. Coulter, Phys. Rev. D51 (1995) 3913; 

H. Dreiner and H. Pois, Two Loop Super symmetric Renormalization Group 
Equations Including R Parity Violation and Aspects of Unification, Preprint 
ETH-TH-95-30 [|iep-ph/9511444| . 



[22] Again, the list of relevant literature is nearly endless. Let me mention, rather 
randomly, a few recent works. Smoothing the d function thresholds in the non- 
supersymmetric theories (e.g. QCD) is discussed in 

D.V. Shirkov and S.V. Mikhailov, Z. Phys. C63 (1994) 463; D.V. Shirkov, 
On Continuous Mass-Dependent Analysis of DIS Data, Proc. Int. Europhysics 
Conf. on High Energy Physics, Brussels, 1995, Eds. J. Lemonne, C. Vander 
Velde and F. Verbeure (World Scientific, Singapore, 1996), page 141; 
J. Chyla and A. Kataev, Theoretical ambiguities of QCD predictions at the Z 
peak, in the Report of the Working Group on Precision Calculations for the 
Z Resonance, CERN Yellow Report 95-03, Eds. D. Bardin, W. Hollik and G. 
Passarino, page 313 ||hep-ph/9502383|[ . 



see also references therein. In the supersymmetric context various smoothing 
procedures are considered in the following papers: 
A. Faraggi and B. Grinstein, Nucl. Phys. B422 (1994) 3; 
M. Bastero-Gil and J. Perez-Mercader, Nucl. Phys. B450 (1995) 21; 
L. Clavelli and P.W. Coulter, SUSY Unification with Smooth Threshold Behav- 
ior, Preprint UAHEP-954 ||hep-ph/950726T|1 . 



[23] S. Dimopoulos and H. Georgi, Nucl. Phys. B193 (1981) 150; 
N. Sakai, Z. Phys. Cll (1981) 153; 

For a review see e.g. R.N. Mohapatra, Unification and Supersymmetry, 2-nd 
Edition, (Springer- Verlag, Berlin, 1992), Chapter 13. 

[24] Perhaps, some elements can be maintained, in a certain limited sense. The hope 
is not yet lost completely. 

[25] D.R.T. Jones and L. Mezincescu, Phys. Lett. B136 (1984) 242; 
A. Parkes and P. West, Phys. Lett. B138 (1984) 99; 
P. West, Phys. Lett. B137 (1984) 371; 

D.R.T. Jones and L. Mezincescu, Phys. Lett. B138 (1984) 293; 

J.E. Bjorkman and D.R.T. Jones, Nucl. Phys. B259 (1985) 533; 

I. Jack, D.R.T. Jones, and C.G. North, N=l Supersymmetry and the three 



24 



loop anomalous dimensions for the chiral superfield, Preprint LTH 368 |[hep- 



ph/9603386 



P.M. Fereira, I. Jack, and D.R.T. Jones, The Three Loop SSM (3 Functions, 
Preprint LTH 372 ||hep-ph/9605440l . 

D. Kazakov, private communication; 
P.M. Fereira, I. Jack, and D.R.T. Jones, The Three Loop SSM (3 Functions, 
Preprint LTH 372 [|hep-ph/9605440|l . 

[27] N. Sakai and T. Yanagida, Nucl. Phys. B197 (1982) 533; 
S. Weinberg, Phys. Rev. D26 (1982) 287; 

M. Vysotsky, Phys. Lett. B114 (1982) 125; T.M. Aliev and M. Vysotsky, Phys. 
Lett. B120 (1983) 119. 

[28] J. Hisano, H. Murayama, and T. Yanagida, Nucl. Phys. B402 (1993) 46. 

[29] B. Blok, Cai-Dian Lu, and Da-Xin Zhang, Proton Lifetime, Yukawa Couplings 
and Dynamical SUSY Breaking m SU(5) GUT, Preprint TECHNION-PH-96-5, 
|hep-ph/96023iq] . 

[30] V. Lucas and S. Raby, GUT Scale Threshold Corrections in a Complete Su- 
persymmetric SO(10) Model: a s (Mz) vs. Proton Lifetime, Preprint OHSTPY- 
HEP-T-95-028 ||hep-ph/ 960 1301] . 



25 



